Preprint typeset in JHEP style - HYPER VERSION 



A Brachistochrone Approach to Reconstruct the 
Inflaton Potential 



Daniel Wohns^*, Jiajun Xu^^^and S.-H. Henry Tye^-f 

-'^Laboratory of Elementary Particle Physics, Cornell University, 

Ithaca, NY 14850, USA 
^Department of Physics, University of Wisconsin-Madison, 

Madison, WI 53705, USA 



Abstract: We propose a new way to implement an inflationary prior to a cosmological 
dataset that incorporates the inflationary observables at arbitrary order. This approach 
employs an exponential form for the Hubble parameter H{(/)) without taking the slow-roll 
approximation. At lowest non-trivial order, this H{(j)) has the unique property that it is 
the solution to the brachistochrone problem for inflation. 



* dfw9@cornell.edu 
^ jxu47@wisc.edu 
*sht5@cornell.edu 



Contents 



1. Introduction 



1 



2. The Minimal E- folds in Canonical Single Field Inflation 



2 



3. The Power Spectrum 



5 



4. Data Constraints 



7 



5. Conclusions and Remarks 



11 



1. Introduction 

In inflationary universe, we can use the evolution of the inflaton as a clock. For a single 
scalar field with canonical kinetic term, the inflationary Hubble parameter can be expressed 
as a function of (p only: 



Since we need an almost flat inflaton potential for enough e-folds of inflation, we may 
choose to start with a truncated expression [1], for some M, 



where <j) is in Planck units and 6j's are dimensionless. The coefficients bi are expected 
to be small, so the truncation to M + 1 terms, with a relatively small M, should yield 
a very good approximation, especially for the region where CMB and other cosmological 
data are available to constrain them. Most efforts to implement an inflationary prior to a 
cosmological data set along this line use H{<f)) of this form [2, 3, 4]. An alternative approach 
uses the inflationary flow formalism to reconstruct models without assuming slow roll [5]. 

One might worry that different parameterizations would lead to different constraints 
on the inflationary parameters. It was shown in [7] that third-order Taylor expansions 
of H{(l)) and H'^{(j)) with a prior of sufficient e-folds yield the same constraints on the 
observable window of inflation. 

The conditions for a reasonable trial H{(j)) are: ffist, it is close to what data indicate, so 
the Taylor expansion needs only a few terms; second, the expression is easy to manipulate 
when calculating observables. Third, it will be nice if the leading term has a simple physical 
interpretation. In the above case, the leading term is simply a flat potential yielding an 
unlimited number of e-folds. -ff(</>) of Eq.(1.2) is just a perturbed version of a flat inflaton 
potential. 





H{^) = Ho{l + bi^ + + • • • + bM4>^^/Ml) 
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We like to ask if there is another expression for H((p) which has these three properties. 
Furthermore, it will be really useful if one does not have to take the slow- roll approxima- 
tion. Here we claim that there is such an expression, which in addition has a conceptual 
underpinning behind it. Consider 

= i/oexp(ai(/> + a20V2 + a3'AV3! + --- + aA/0^V^!) , (1-3) 
here H^^\(f)) = Ho corresponds to a flat potential. What is particularly interesting about 
this exponential form is that, when it is truncated to the linear term in the exponent, i.e., 

= Hoexp{ai())) = Ho exp{^/^ (/)) (1.4) 

it is simply the solution to the brachistochrone problem for inflation without taking the 
slow-roll approximation. That is, away from the slow-roll approximation, H^^\(p) yields 
the minimum number of e-folds (fastest path) for a fixed drop in H over a fixed field range. 
So, any deviation from this path will yield more e-folds. In this sense, it is the opposite of 
the fiat potential. 

The exponential form of H((j)) was also mentioned in an early paper [6] by Andrew 
Liddle, as a way to construct arbitrary inflaton potential. Here, we motivate the form 
Eq.(1.3) from entirely new perspective, i.e. the brachistochrone problem for infiation and 
how to deviate away from the brachistochrone solution. In this sense, we have given the 
exponential of H{(p) an interesting physical interpretation. 

Furthermore, it is quite amazing that, for small e, H^^\(f)) alone (i.e., without higher 
terms in Eq.(1.3)) can yield an infiationary scenario within experimental bounds, implying 
a low (small M, say M > 2) H^^'^^cj)) will do very well already. Models that are close to 
the brachistochrone solution typically require a spectral index that is within a few percent 
of unity. Lastly, the relations between the slow-roll parameters and the coefficients aj in 
H^^\(f>) of Eq.(1.3) are very simple. 

2. The Minimal E-folds in Canonical Single Field Inflation 

The dynamics of a single canonical scalar field during infiation can be described by the 
Hamilton- Jacobi formalism [8] , where we write every quantity as a function of the scalar 
field (j), i.e., we choose cj) as the "clock". We will use the notation H^j, = dH/dcj), H^^ = 
(fH/d(j)'^, etc. 

3H^ = ^4>^ + V{^) , (2.1) 
= 2H^{(P) . (2.2) 

We set 'rn~^ = 8ttG = 1 throughout the paper. Here we do not consider violation of the 
Null Energy Condition, so H < 0. We choose the convention (/) < so that dH/dcj) > 0. 
The e and r] parameters in this formalism are given by 



The number of e-folds is given by 

Consider the case with fixed initial Hubble scale Hq = H{(j) = 0) and final Hubble scale 
Hf = H{—A4>) over a fixed finite range A(p of (f). We can ask what kind of function H{(j)) 
gives the minimum number of e-folds between these two fixed boundary points. Starting 
with Eq.(2.5), one can write down the Euler-Lagrangian equation, 



H\ d d f H 



dH\H^I dd) dH^ \ Ha 

which gives 



(2.6) 



H \H 

Comparing with the definition of r] in Eq.(2.4), the Euler-Lagrangian equation is equivalent 
to saying that r] = for the trajectory with minimum e-folds. 
The solution to Eq.(2.7) is given by 

H{^) = Hoexp{ai<l)) (2.8) 

where the integration constant oi is related to the e parameter as 



2al 



For inflation, we require 



0<e<l ^ 0<ai<^. (2.9) 



The lower bound ai > arises because dH/dcp > 0. 

From H{(j)) above, we can reconstruct the potential through Eq.(2.1) and Eq.(2.2), 

V {(/)) = {3 - e) exp(y2~e(l)^ . (2.10) 
In the limit e ^ 1, the potential becomes linear 

Vs{(l))^3HUl + V2eA , (2.11) 



which agrees with previous analysis in the slow-roll scenario with large damping [9] ^. 
One can now calculate the value of the minimum number of e-folds 

iV- - ^ / - (2 12) 

- 2 J "^^H^ic^)- 2ln{Ho/Hf) ' ^''^'^ 



^Note that, to get the correct power spectrum mdex, we have to keep the 0^ term in Vs{<j))- Using 
Eq.(3.3) one obtains, Us — 1 = dlnPTj/dlnfc = ^(2e + r])/{l ~ e) ~ — 2e — rj ~ + 2iisr where the last 
equality is taken in the slow-roll approximation. 
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Since the minimal e-fold path corresponds to 77 = 0, we can turn on the 77 parameter to 
characterize the deviation away from the minimal e-fold path. With non-zero rj, we need 
to solve 



IT 

The solution can be formally written as 



H 4 



H( 



Hq exp 



We can fix the integration constant a[ by matching H{—A(f)) = Hf. Therefore 



HHf/Ho 



A(j) 



-A(j> 



Now the number of e-folds is 



H 



-d</. 



A<t> 



d(/>2??(02 



(2.13) 



(2.14) 



(2.15) 



(2.16) 



In principle, one can consider an arbitrary function of r/(0), as long as |7/((/')| < 1- Here 
we will illustrate using a constant rj, so that the solution Eq.(2.14) can be simplified, 



H 



Hq exp ya'icl) 
HHf/Ho) 



In 



7] 



In the small r] expansion, we can write 

(A0)2 



+ 



Act) 

HHf/Ho) - |7?(A</>)2 

HHf/Ho) + ivi^^r 

{A(P)S^ 



+ 0(r?*) + 



(2.17) 
(2.18) 

(2.19) 
(2.20) 



2\niHo/ Hf) 38i[HHo/Hf)]^ 

which further verifies that — t- 0, Ng approaches the minimum. 

Instead of the ansatz (1.3), we can also use the following ansatz 

H{(l)) = Ho e""^ (1 + bi<P + b2<p^/2 + ...) (2.21) 

where 

e = 2(ai + 6i)2 (2.22) 

This interpolates between the two limiting starting functional forms for the inflaton poten- 
tial. If we want to increase by one parameter only, we can simply set 61 = 0, 

H{(t>) = Ho e""^^ (1 + b2<p^/2) (2.23) 

Now 

e = 2(ai + 6i)2 = 2a? , r/ = 4(6? - 62) = -462 • (2.24) 

Starting from H^^\(j)), it is unclear which is the best ansatz to generalize it, Eq.(1.3) or 
Eq.(2.21). It will depend on future data. 
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3. The Power Spectrum 



We have seen that the form of H{(j)) we proposed has some theoretical underpinning: 
H^^\(j)) is the form that gives the minimum number of e- folds, and H^'^\(j)) can be inter- 
preted as deforming away from the minimal e-fold trajectory using a constant r] parameter. 
Going beyond H^^\ it is not clear what is the best way to generalize the form of H{<j)). At 
least we have seen two possibilities, Eq.(1.3) and Eq.(2.21). In this section, we will try to 
argue that the form Eq.(1.3) is more favorable in the sense that inflationary observables 
can be more simply expressed in terms of its parameters. We will see that the tensor 
and scalar power spectrum, their spectral index, running of spectral index, etc., can all be 
written entirely in terms of \nH and its (f) derivatives, which makes the exponential form 
of Eq.(1.3) extremely convenient for computing observables. 
Let us first look at the scalar power spectrum 

Pn{k) = ^= rr,2 . „ > (3-1) 



87r2e 16TT^{ln H)l 



<f>=0 



where, without loss of generality, we have chosen (p = to correspond to the moment when 
the wave vector k crosses the horizon. 

Furthermore, we notice that the e and r] parameters can be written in terms of In H 
and its derivatives 

e = 2{lnH)l, = -A{\i^H)^^ . (3.2) 
We also have the derivative 

d ^ d(/> d ^ -2(lng)^ d 
dlnfc dlnA:d(/) 1 - e d(^ ^ ' ' 

We see that In if appears everywhere, which suggests that if we Taylor expand hiH 
instead of i.e., using the exponential form in Eq.(1.3), the resulting din A; derivatives 
will be very easy to compute. We now show explicit computations of n^, dn^/dlnfc, r and 
nt starting with H{(j)) of Eq.(1.3). 



din /c 



A{\nH)^^-A{\uH)l 



<t>=o 



l-2{lnH); 

4a2 - 4af 
1 - 2af 



0=0 



(3.4) 



so the scalar spectral index is 



= 1 + J^f p. 1 + 4a2 - 4a? (3.5) 

To compare with the e and rj expansion, simply replace e = 2af , ry = — 4a2, and we get the 
usual result Ug = 1 — 2e — rj. The advantage of writing H{(j)) in exponential form is that 
by derivatives of In//, one can extract the coefficients Oj's more transparently. 
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Comparing to the usual e and rj expansion, here the order of expansion can be tracked 
by counting the total number of (p derivatives. So af and 02 are of the same order, since 
they measure the two ^ derivative terms. The parameter as is of higher order, since it 
measures the term with one more derivative of (f). 

To further appreciate the advantage of the exponential form, let us compute the run- 
ning of the spectral index. Up to four derivatives, we get 



dns 
dink 



-2{lnH)^ d 
1-e d(t) 



[4(lnif)^^-4(lni?) 



4>=o 



-80103 + 16afa2 
l-2af 



(3.6) 



In a more careful analysis including 0{e) corrections in the power spectrum. 



Pn 



l_(C + l)e-icr? 



n 2 



C = -2 + ln2 + 7 



167r2(hii7)2 

where 7 is the Euler-Mascheroni constant, we get up to four derivatives, 

-4(lni?), 



Us 



1 + ^?^ + 



2al 



1 



4(C + l){\iiH)^{\nH)^^ + 2C{\nH) 



40,2 - 4of 16(C + l)of 02 - SCoiOa 



2o2 



2o2 



(3.7) 



and up to six derivatives 



dus _ -2{\ilH)^ d 
din A; ~ 1 - 2(lni7)2 d^ 

+16(C + l)(lniJ)2(lnii-)^^ - 8C(lni7)^(lniJ) 



160^02 — 80103 
l-2of 



+ 



1 - 2of L 



64oto2 - (32C + 48)a?03 



-(64C + 96)o?oi + 16Co?04 + 16Cai0203 

Similarly, for the tensor power spectrum 

2/^2 



(3.8) 



Pt = 



[1 - (C + 1)6] 



we get the tensor to scalar ratio 



and the tensor spectral index 

-2{\nH), 



nt 



l-2(Ini/)2 
-4of + 16(C + l)ofa2 



32(lni7)^(l + Cr?) 
32o?(l-4Co2) , 



[2(lni7)^-8(C + l)(ln/f)^(ln/i' 



(3.9) 



1 



2o2 



(3.10) 
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Our approach can compute higher derivatives d^Us/ {din k)^ (n > 3) very efficiently, 
as all the derivatives of InH are tied to its Taylor expansion coefficients directly. The 
computational effort does not increase much with taking more derivatives, and new pa- 
rameters 04 , as , . . . , enter in a systematic way. The total number of (p derivatives in 
each term can be used as order counting in the expansion. Whereas in the usual approach, 
the higher derivatives of the expansion parameters e, r/, ^ are tedious to work out. Clearly, 
the exponential form of H{(j)) has its advantage in computing power spectrum observables. 

4. Data Constraints 

We use Markov Chain Monte Carlo methods to fit WMAP7 data, using modified ver- 
sions of CAMB [11] and COSMOMC [12]. The independent parameters we choose are 
{01,02,03,04}. All other cosmological parameters are fixed at the central values in the 
ACDM fit of the WMAP7 data in Ref. [13]. We parametrize the scalar and tensor primor- 
dial power spectrum by 

(T \ ns-H-|(dns/dlnfc)ln(fc/feo) 

Pr = rA. (-) (4.2) 

The relation between {77,^, dn^/dlnA;, r} and {01,02,03,04} are worked out in Sec. 3. The 
scalar amplitude is fixed at As = 2.23 x 10~^ by adjusting the parameter Hq, i.e. given oi 
and As, we set Ho/mpi through the relation 

As = ^0 . 

^ levr^ofmpj 

Since we fix ^4^, oi is correlated with H/mp\. For inflation, we require < oi < ^/2/2. In 
our data analysis, we put a lower bound on oi by hand, i.e. we searched the parameter 
space with 10~^ ^ < a/2/2. 

The theoretical priors we impose are 

• Perturbative expansion of \n{H): We require that higher-order terms in the Taylor 
expansion of hi{H) are smaller than leading-order terms, that is 

ai4> > a2(l>^ , aicj) > 03^^ , ... (4.3) 

• Monotonic function for H{(f)): During inflation, (f) cannot change sign, which implies 
that is a monotonic function. In our convention, we require ff^ > 0. For the 
cubic truncation of In(ff), 04 = 0, H^f, > is equivalent to the following constraints 
on oi, 02 and 03 

03 > , 03 - 2ai03 < . (4.4) 

If 04 7^ but small, the above constraints will be relaxed, but we do not expect 
qualitative change of the parameter space. As we will see, data currently are not 
sensitive to 04. In our analysis, we impose the same constraint (4.4) even if 04 7^ 0. 
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Figure 1: One-dimensional marginalized constraints of parameters from WMAP7. {01,02,03,04} 
are free parameters and the rest are derived using equations in Sec. 3. Solid lines show the full 
marginalized probability density on a linear scale. Dashed lines show the mean likelihood of samples 
on a linear scale. 



The end of inflation: We assume there are two ways to end inflation in our approach: 
(1) Inflation ends when e grows and becomes larger than unity. (2) Although e does 
not grow larger than one during inflation, when the inflaton travels for A.(f) ~ 1, 
higher-order terms in H{(j)) becomes important; they may change the form of the 
potential to shut off inflation. 



One can estimate how much the inflaton field travels before e grows to one. This 
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gives 



Ac/) 




al + -v/2a3 — 2aia3 



(4.5) 



Note that with > and oi < \/2/2, we have a| + V^a-^ — laia^ > 0, so Acf) is 
always well defined. 



Therefore, inflation ends with e = 1 if 



a2 + J a\ + \/2a3 - 2aia3 



< 1 . 



(4.6) 



Otherwise, inflation ends when A.(j) = — 1. 

• Number of e-folds: We require that after the fiducial scale fco = 0.002 Mpc~^ has left 
the horizon, inflation will last for A^e e-folds. Since we are only interested in large 
scale observations, we constrain Nf, in a small window 



By imposing the lower bound on N(,, we discard models which do not yield enough 
e-folds to solve the horizon problem. The upper bound on N(. makes sure that our 
parameters {01,02,03,04}, which are defined at = 0, correspond to scales within 
the observable window of WMAP data. 

Fig.l and Fig. 2 show constraint on our parameters {oi, 02, 03, 04} from WMAP7 data. 
We have run eight MCMC chains and they converge with the Gelman and Rubin statistics 
R — 1 < 0.01. Since we only use the data analysis for illustration, we have not run the 
MCMC for longer time to achieve better convergence. 

We see that data constrain the parameter oi to be around oi ~ 0.006, corresponding 
to e ~ 10~^. This is consistent, since we are essentially dealing with small field inflation 
with Ai;^ < 1. The total number of e-folds are N(, ~ l/-v/e ~ 10^. 

The lower bound on oi may come as a surprise at first sight. One may think that 
given a reasonable value of oi, we can always change the value of H/nipi to fit the scalar 
amplitude, so there should not be a lower bound on oi as shown in Fig.l. In fact, the lower 
bound comes from our theoretical prior > 0, i.e. Eq.(4.4), which leads to oi > 02/(203) 
for given 02 and 03. One may increase 03 to allow smaller value of oi, but on the other 
hand, our perturbative truncation requires that a^cj)^ < aicp, so 03 cannot be increased 
without limit. In summary, the constraint on oi mainly comes from the set of theoretical 
priors we applied in the data analysis, in particular our choice A(j) = — 1. 



A^pivot - 3 < Ae < Apivot + 3 



(4.7) 



while Apivot is correlated with H/nipi through 




(4.8) 
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As e ~ 10"'^, its contribution to Ug is negligible. So is determined mostly by 02- That 
is why 02 becomes the best constrained parameter among the four. Currently, WMAP7 
data do not constrain 03 and 04 very well, which suggests that the quadratic truncation 
i?((/>) in Eq.(1.3) already provides a very good fit of data. 

Our analysis gives a slightly blue spectral index with a small negative running param- 
eter. This is consistent with the WMAP7 analysis [14]. When we allow the spectral index 
to run, the central value of Us shifts towards the blue-tilted side. However, one should 
not take seriously our result of a nonzero dng/dln/c. This is most likely an artifact due 
to fixing other cosmological parameters. For example, Ref.[4] has shown that the value 
of total matter density is strongly correlated with the running of spectral index. We 
expect that when marginalizing over all ACDM parameters, the running of spectral index 
will likely become insignificant. 

At fixed order, our theoretical priors impose an upper bound on jn^ — 1|. At second 
order this bound is, using Eqs. (3.5) and (4.3): 



4ai(A(/))-i + 4a? _ 2^/2e(A(/.)-i + 2e 

At third order this bound is: 



ns-l\< — ^ = — — (4.9) 

' l-2a? 1-e ^ ^ 



I , , ^ 2V2^(A,/.)-i + (2 - 4C(A</>)-^)e + 4^/2(C + l)e3/2(A<^)-i 

\ns - 1| < (4.10) 

where C ~ —0.7. This bound gets progressively weaker at higher order. For lS.cj) = 1 and 
e <C 1 at low order, 

\ns - 1| < 2^2^. (4.11) 

For e « 10~^, Eq. (4.11) implies \ns — 1| ^ 0.03. A small perturbation of the brachis- 
tochrone solution (1.4) with Ac/) = 1 is only a good fit to data when Eq. (4.11) is satisfied. 
Note that in deriving this bound, we assumed that the leading order term in (1.3)is domi- 
nant and that the form (1.3) is valid throughout the final sixty e-folds of inflation. Relaxing 
one or both of these assumptions would eliminate our theoretical prior and allow the spec- 
tral index to deviate further from unity, but render the brachistochrone interpretation that 
originally motivated our approach inapplicable. 

The constraints on r are driven primarily by the constraints on oi. That means the 
reported bounds on tensor to scalar ratio is primarily from model, i.e. the inflationary 
prior, not from data. The bound on r from data is much weaker according to the WMAP 
analysis [14]. In our analysis, we focus on the set of models with A.(f) < 1, which will limit 
the level of tensor mode according to the Lyth Bound A(j)/ ANe ~ \/r/8. Our method can 
also be applied to cases with Acp > 1, but then one have to worry about the truncation 
of higher order terms in the expansion of In H. We will leave such an analysis for future 
work. 

It will be very interesting to include more data sets in our analysis, such as ACBAR, 
Large Scale Structure and PLANCK data. With a wider span of scales, we may be able 
to achieve better constraints on 03 and 04. In addition, PLANCK will provide us with a 
tighter bound on tensor mode, which will help to further constrain oi. Furthermore, given 
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our exponential form of H{(j)) and the polynomial form in Ref.[l, 4], it will be interesting 
to see which functional form is favored by data, and to see if single-field models can be 
distinguished experimentally from more general models [15, 16]. Since PLANCK and other 
data sets will become available soon, we leave the full data analysis and Bayesian model 
comparison for future work. 



5. Conclusions and Remarks 

In this paper, we have proposed a new functional form of H{(j)) as an effective way to 
implement inflationary priors in cosmological data analysis. This new form of H((j)) entails 
writing H in the exponential form. At the lowest order, with the exponent linear in it is 
the solution to the brachistochrone problem for inflation, which corresponds to the minimal 
number of e-folds for a flxed drop in H over a flxed fleld range. Higher-order terms can be 
included to deviate the inflaton path away from the brachistochrone solution. 

In addition to the theoretical underpinning, the exponential form of H{(j)) also provides 
an efficient way to compute power spectrum observables. If one Taylor expands the function 
ln{H{(j))), the expansion coefficients are natural parameters to express the observables, such 
as Hs, dng/dlnk, r and nj. Higher-order derivatives of the power spectrum entails higher- 
order expansion coefficients of ln(i7(i;^)), and the computation is more straightforward than 
the usual e and r] parametrization. 

We have also performed MCMC analysis to illustrate how observation data can con- 
strain the expansion coefficients of ln{H{(p)). The WMAP7 data constrains the two leading 
coefficients very well, but is not quite sensitive to the coefficients of the cubic and quar- 
tic terms. We expect including more data on different scales will improve the constraints. 
With the upcoming precision measurements from Cosmic Microwave Background and Large 
Scale Structure, we hope our proposal will offer an efficient way to reconstruct the inflaton 
potential from data. 

In the actual implementation of an inflationary prior to a cosmological dataset, one 
does not need to use the same H((p) for the whole range of sixty e-folds of inflation. The 
actual inflaton potential can be more complicated; in fact, it can be multi-fleld. In these 
more general situations, one can use a set of piecewise exponential segments of the form 
(1.3) instead. In this case, the constraint of sixty or more e-folds imposed in the above 
analysis may be relaxed. 
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Figure 2: Two-dimensional marginalized distributions from WMAP7 data. The contours represent 
68% and 95% confidence levels. The color represents mean likelihood of the samples with red being 
high likelihood, followed by yellow. 



- 13 - 



